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ABSTRACT 

We report on the mass-to-light ratio determination based on a newly selected 
binary galaxy sample, which includes a large number of pairs whose separations 
exceed a few hundred kpc. The probability distributions of the projected 
separation and the velocity difference have been calculated considering the 
contamination of optical pairs, and the mass-to-light ratio has been determined 
based on the maximum likelihood method. The best estimate of M/L in the B 
band for 57 pairs is found to be 28 ~ 36 depending on the orbital parameters 
and the distribution of optical pairs (solar unit, Hq = 50 km s~^ Mpc~^). 
The best estimate of M/L for 30 pure spiral pairs is found to be 12 ~ 16. 
These results are relatively smaller than those obtained in previous studies, but 
consistent with each other within the errors. Although the number of pairs with 
large separation is significantly increased compared to previous samples, M/L 
does not show any tendency of increase, but found to be almost independent of 
the separation of pairs beyond 100 kpc. The constancy of M/L beyond 100 kpc 
may indicate that the typical halo size of spiral galaxies is less than ~ 100 kpc. 
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1. Introduction 



The mass of a galaxy is a fundamental quantity in understanding its dynamics and 
structure. Mass distribution in galaxies has been extensively studied with optical and HI 
rotation curves. Several studies revealed that spirals galaxies have flat rotation curves even 
at the observed outermost points, indicating the existence of extended dark halos (e.g., 
Sancisi & van Albada 1987) . The extent and total mass of dark halos are, however, not 
understood well and yet to be studied in detail. For further investigation of extended halos. 
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different approaches are required to trace the mass distribution beyond the HI disk, where 
rotation curves cannot be measured. 

Binary galaxies are useful for the determination of the total mass or mass-to-light 
ratio (M/ L) of galaxies, like stellar masses are measured from the motion of binary stars. 
Unlike stellar binaries, however, the total mass of individual binary cannot be directly 
determined because of the long orbital periods. Instead, statistical treatment is necessary 
to obtain the average mass or M/L of the sample galaxies. Many efforts were made to 
determine the M/L ratio of binary galaxies statistically (e.g., Page 1952; Karachentsev 
1974; Turner 1976a, b; Peterson 1979a, b; White 1981; van Moorsel 1987; Schweizer 1987; 
Chengalur et al.l993; Soares 1996). In binary galaxy studies, a careful selection of binary 
galaxies is very important, since the biases in selecting pairs should be corrected for to 
determine the mass or M/L ratio. Turner (1976a) proposed well-defined selection criteria 
based only on the positions and magnitudes of galaxies. Later investigations (e.g., Peterson 
1979a) also made use of similar selection criteria independent of radial velocities, which 
are so-called 'velocity-blind' pairs. According to such velocity-blind selection criteria, 
two galaxies are regarded as a pair if they have no close companion compared to their 
projected separation. This 'velocity-blind' selection criterion is simple and convenient for 
pair selection, but its problem is that the criterion could introduce strong bias toward pairs 
with small separations; for pairs with wider separations, company galaxies are searched for 
in a larger region, leading to exclusion of widely-separated pairs with higher probability. In 
fact, the average separations of selected pairs in these studies were 50 kpc ~ 100 kpc (see 
Peterson 1979b). Since the dark halos could extend beyond this range, it is important to 
study binary galaxies further based on pairs with wider separations. 

The other major problem of the velocity-blind sample is that the sample suffers from 
the contamination of 'optical pairs', which consists of two isolated galaxies projected close 
by chance. In order to reduce the contamination of optical pairs, it is better to select binary 
galaxies based not only on the positions but also on the radial velocities. Fortunately, the 
number of radial velocity observations is rapidly increasing thanks to recent large-scale 
redshift surveys. Moreover, the observational uncertainty has been significantly reduced 
due to the recent development of observational instruments, which enables us to estimate 
M/L with better accuracy compared to previous studies. Therefore, it is interesting to 
study binary galaxies again by utilizing such huge data. 

For these reasons, in this paper we study the mass-to-light ratio of binary galaxies 
by making use of the database. The plan of this paper is as follows. In section 2, we will 
describe how to select widely-separated pairs effectively, while reducing the contamination 
of optical pairs. The selection criteria, and the basic data for selected pairs will be presented 



-3- 



in section 2. In section 3 we will perform maximum-likelihood analysis based on the orbital 
models of binary galaxies, and determine M/L. We will also consider M / L dependence on 
galaxies' type. Discussions on the dark halo extent will be given in section 4. 



2. Selection of Pairs 

2.1. Basic Idea, Selection Criteria, and Sample 

The observable quantities for orbital motion of pairs are the projected separation rp, 
and the radial velocity differences Vp. A set of and fp can be used to estimate the total 
mass of a pair through an estimator of mass, for example, r^v'^/G. However, the mass of 
galaxies varies by about 3 order of magnitude from dwarf galaxies to giant ellipticals. A 
better quantity which represents the mass of galaxies is the mass-to-light ratio, M/L. The 
mass-to-light ratio of galaxies is expected to vary much less than the mass itself, and hence 
we focus on the mass-to-light ration of pairs in this paper. 

What can be obtained through binary galaxy analysis is the total mass to total hght 
ratio of pairs which is written as (Mi -|- M2) / {Li + L2), but in the rest of this paper we denote 
this ratio as M/L for simplicity. Note that if Mi/Li — M2/-L2, the total mass-to-light 
ratio M/L is equal to Mi/Li and M2/L2. For convenience in M/L estimate, we define the 
luminosity-corrected separation and the luminosity-corrected velocity difference as 

i?p ^ rp/Ll/^ (1) 

T/p ^ W\/L'l\ (2) 

A combination of i?p and can give an estimator of the mass-to-light ratio of pairs. This 
estimator, which we call the projected mass-to-light ratio, is defined as 

{M/L), ^ ^ = (3) 

If a bound pair of galaxies are separated so widely that they can be approximated as two 
point masses, the law of energy conservation gives that 



2 
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< M/L, (4) 



because the total energy for a bound pair is always negative. A combination of equation (3) 
with inequality (4) gives 

(M/L)p < 2(M/L). (5) 
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If all pairs are bound and have the same M/L, the binary populations in the R^-V^ phase 
space lie below the envelope which corresponds to 2(M/L). In practice, the number of 
pairs is limited and insufficient to see the true envelope corresponding to 2[M/L) in the 
Rp-Vp space. Detailed calculations of probability distribution show that pairs are likely to 
concentrate to small Vp, and thus smaller {M/L)p due to projection effect(e.g., Noerdlinger 
1975). We will discuss this in later section by calculating the probability distribution based 
on the Monte-Carlo simulation. In any case, the pair distribution in the Rp-Vp phase space 
can be used for testing whether or not bound pairs are efficiently selected: while bound 
pairs are hkely to have small Vp, optical pairs could have extremely high Vp and {M/L)p. 

Here we describe the selection criteria for pairs. For convenience, we define the total 
luminosity of a pair normalized with lO^^L© in the B band as Lio = {Li + L2)/(1O^°L0). 
Note that this luminosity roughly corresponds to that of the Milky Way Galaxy. We 
re-define Rp and Vp normalizing with luminosity Lio as, 

= 7173 (kpc), (6) 

and 

^p^Mlknis-^). (7) 

As the first step of pair selection, we have to select a pair of galaxies that are relatively 
close to each other both in the sky plane and in the redshift space so that they are likely to 
be bound. Since the average separations in previous studies were around 100 kpc, and since 
we are interested in widely-separated pairs, we set the maximum projected separation of a 
pair in the sky plane as 

1: Rp< 400 kpc. 

The maximum velocity difference must also be large enough to include pairs orbiting 
around each other at high velocity. Since galaxies as luminous as 1O^°L0 have rotation 
velocity of ~ 200 km s~^, the velocity difference of a pair could be as high as a few hundred 
km s~^. Hence, we set the maximum velocity difference of pairs as 

2 : Vp < 400 km s'K 

Note that the radial velocity difference is usually quite small compared to the true velocity 
difference due to projection effect (Noerdlinger 1975). Thus, physical pairs are unlikely to 
have larger radial velocity differences than the maximum value given above. 

The observational data set cannot be complete to faint galaxies, and hence we limit 
the application of our analysis to sufficiently bright galaxies. Since nearby galaxies are 
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cataloged almost completely down to 15.5 magnitude (e.g, de Vaucouleurs et al. 1991), we 
set a criterion for the B band magnitude as 

3 : 7711,1712 < 15.0 mag, 

and also set an upper limit for the total magnitude of a pair in the B band, mi+2, as 

4 : mi+2 < 13.5 mag 

In order for selected pairs to be likely to be bound, pair galaxies must be isolated 
well. We regarded two galaxies as a pair if all of its companion galaxies brighter than 
mi+2 + 2.0 mag satisfy both 

5 : ^ > a 400 kpc, 

and 

6 : ^ > 6 400 kms-^ 

Here ri,and Vi are the projected separation and the radial velocity difference of ith 
companion galaxy with respect to the luminosity center of the pair. Note that we set the 
lower hmit of the total magnitude mi+2 to be 13.5 mag (criterion 4). The faintest galaxies 
that should be considered is, hence, at 15.5 mag, to which magnitude galaxies are cataloged 
almost completely. 

Parameters a and b determine the volume for companion search, and so determine the 
degree of isolation. Note that the volume depends only on the total luminosity of a pair, 
Lio, but independent of the separation Rp of a pair. Therefore, as far as pairs with the 

same luminosity are concerned, companions are searched for in the same volume, and thus 
the criterion does not introduce bias toward pairs with small separations. We set b = 1.5 
throughout this paper, but tested three values of a (1.5, 2.0, and 2.5) to seek a value for 
effective selection of bound pairs. 

We applied criteria described above to the sample of galaxies that we compiled for this 
study using NED (NASA Extra-galactic Database). The sample consists of bright nearby 
galaxies with redshift less than 4,500 km s^^ . The upper limit for redshift is introduced 
because the number of bight pairs which satisfy the criteria 3 and 4 becomes small at large 
redshift. The data for positions, heliocentric velocities, and the B band magnitudes were 
mainly taken from NED, and supphed with RC3 catalog (de Vaucouleurs et al.l991). The 
distances to the pairs are obtained using the redshift of the luminosity center {Hq — 50 
km s"^ is assumed). In order to avoid the error in distance due to the local deviation 
from Hubble flow, galaxies with redshift smaller than 1,000 km s~^ are excluded from the 
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sample. Galaxies in clusters and close to clusters may deviate from the Hubble flow even 
beyond the rcdshift of 1000 km , but this effect is expected to be small because the pairs 
selected with criteria 5 and 6 are likely to be field binary galaxies. Galaxies with b < 20° 
are also excluded from the sample since objects at low galactic latitude may be significantly 
obscured by galactic extinction The sample we compiled consists of 6475 galaxies with 
magnitude brighter than 15.5 mag and rcdshift between 1000 km s~^ and 4500 km s~^. The 
uncertainty in the magnitude is typically 0.2 mag., which leads to the uncertainty in L of 
20 %. The corrections for intrinsic absorption and galactic extinction were made according 
to de Vaucouleurs et al.(1991). 

The sample is, of course, incomplete in terms of redshift because rcdshift measurements 
were not made for all galaxies. This incompleteness leads to possible mis-identification of 
pairs, if the criteria described above are applied only to a sample of redshift-know galaxies. 
To correct for the effect of redshift incompleteness, primary binary candidates are at first 
searched in the sample of redshift- known galaxies, and then a redshift-blind search was 
performed for the primary binary candidates. If there is any redshift-unknown companion 
which is brighter than mi^2 + 2.0 mag and is so close to the pair that the criterion 5 is 
violated, the pair was rejected from the binary candidates. About 30 % of pairs in the 
primary binary candidates were rejected through this procedure. 

The sample of binary galaxies after the correction for the redshift incompleteness 
still contains some pairs that are not appropriate for this study. For instance, the basic 
data for the analysis, such as Vp, mi, and m2 could be quite uncertain for some pairs. 
In particular, the uncertainty in the radial velocity is crucial for M/L determination, as 
the M/L estimator depends on V^. Therefore, if redshift uncertainty is not reported for 
any of the two galaxies of pair, the pair is excluded. This process reduced the number of 
pairs by 7%. If the magnitude uncertainty and the absorption-corrected magnitude are not 
available, the pair is also excluded, and in this process 7% of primary binary candidates 
were rejected. Moreover, a galaxy could appear in the binary sample twice or more with 
different partner. This can happen if one of pair galaxies is a bright galaxy like the cD 
galaxy and it has several companion galaxies around it. In this case, however, these galaxies 
should be regarded as cluster or group rather than binary galaxies. Therefore, wc also 
excluded possible clusters or groups of galaxies that appear in the binary sample twice or 
more. We found only two possible groups in the primary binary candidates. 
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2.2. Results 

Figures 1 show the distribution of thus selected binary galaxies in the Rp-Vp phase 
space. Two cases for the isolation parameter, a=1.5 and 2.5, are shown. The number of 
selected pairs is 109 and 57, respectively. In the case of a=1.5, the pairs in the -Rp-Vp space 
shows only weak concentration toward small Vp, and a large number of galaxies have high 
(M/L)p exceeding a few hundred Mq/Lq. Even if their true M/L is a few hundred, it is 
unlikely that so many galaxies appear to have so large {M/L)p in the projected phase space, 
as {M/L)p is expected to be significantly smaller than the true M/L due to projection 
effect (Noerdlinger 1975; see also Section 3 of the present paper). This indicates that they 
are probably optical pairs, and that the degree of isolation is not strong enough to select 
bound pairs effectively. On the other hand, for a = 2.5, the concentration of pairs to Vp=0 
is much more clear than for a = 1.5. Most of 57 pairs are distributed below (M/L)p of 20 in 
solar unit, and there are only few galaxies that have high (M/L)p. This correlation between 
Rp and Vp are naturally explained if the separation of pairs are larger than the extent of 
halos so that the pairs can be approximated as point masses (but note that even in case of 
extended hale such an envelope would appear in the projected phase-space like point-mass 
cases; see Soares 1990). In the rest of this paper, we use the binary galaxies sample selected 
with a = 1.5 and 2.5 for M/L determination. We call the sample selected with a = 2.5 as 
sample I, and the one selected with a = 1.5 as sample II. Table 1 summarizes the basic data 
for 57 pairs in sample I. 

3. M/L determination 

In this section, we estimate the M/L ratio of the sample pairs selected above. We 
construct the orbital models for physical pairs, and calculate the probability distribution of 
pairs in the Rp-Vp phase space considering the contamination of optical pairs. Then, we 
compare the models with the observational data, and determine the M/L ratio based-on 
maximum-likelihood analysis. 

3.1. Distribution of Bound Pairs 

First we construct models for orbital populations of binaries. For simplicity, binary 
galaxies are treated as point masses in the following analysis. As can be seen in figure lb, 
pairs show strong concentration toward small Vp in the Rp-Vp space, which is just expected 
from the point-mass assumption. Further tests for validity of the assumption will be made 
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in the next section. 

An ensemble of well-mixed binary population satisfy the Jeans equation (Binney and 
Tremaine 1987), _ _ 

dr r 

Here u denotes the separation distribution of pairs, and /? is the anisotropy parameter 
defined as 

where vq and denote each component of velocity ellipsoids. Note that P — — oo for 
circular orbits, /3 = for isotropic orbits, and /3 = 1 for radial orbits. We may rewrite 
equation (8) by normalizing with luminosity as 

djiyW) 2v0V? _ GM 
dR ^ R R-'L ' ^ ' 

where R — r/V-/^ and = vljl?!'^. Note that they are true separation and velocity, but 
not projected ones. 

For the separation distribution z/, we also assume a power law with an inner cutoff 
radius rmin as, 

v{R) (X R-^ for R > R^^^. (11) 

We introduced the cutoff radius because galaxies have finite sizes, and pairs that are too 
close are not hkely to exist. The model used here is, therefore, not exactly a scale-free 
model (cf. White 1981). 

In order to model the distribution of pairs, one should choose suitable values for 
parameters 7, and Rmm- The parameters related to the separation distribution are 
obtained directly from observed separation distribution, because the probability distribution 
for projection effect can be written analytically as 

P[^pl^] - ^R^R2^R2y/2 (for ^P<^)- (12) 

We compared the separation distribution of observed and model pairs, and obtained 
the best-fit values 7 = 2.6 and i?min = 10 kpc. In the rest of this paper, we adopt these 
best-fit values for 7 and -Rmin, but we note the results are not sensitive to changes in 
the assumed values. Once the separation distribution is obtained, the distribution of 
the velocity difference is obtained by solving the Jeans equation [eq.(8)]. Then, one can 
calculate the probabihty distribution of pairs in the Rp-Vp phase space, Phin[Rp: Vp\{M/L)], 
by taking the projection effect into consideration. 
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3.2. Selection Bicis and Contamination of Optical Pair 



In addition to the orbital models for true pairs, here we consider the selection effect 
for pairs and the contamination of optical pairs. As described in the previous section, 
the isolation criteria are independent of the separation or velocity difference of pairs, and 
hence, the sample is free from biases both for the separation in the sky plane, and for 
the separation along the line of sight. The isolation criteria, however, may cause possible 
exclusion of true pairs due to chance projection of another companion galaxy that are not 
physically related to the pair. According to the isolation criteria, the maximum velocity 
difference of a pair is ~ 400 km s^^ for galaxies with L ~ 10^°Lq. Therefore, a system of 
a true pair plus any foreground or background galaxy at distant within 8 Mpc from the 
pair cannot be regarded as a true pair because the criterion 6 is violated {Hq = 50 km 
s"^ Mpc~^ assumed). Unfortunately, there is no way to determine whether the observed 
velocity difference of two galaxies is due to Hubble flow or due to binary orbital motion, and 
so this kind of exclusion of true pairs is unavoidable. Furthermore, two galaxies which are 
separated well along the hne of sight and are not physically associated could be regarded 
as a pair because of misidentification of the redshift as binary orbital motion. For these 
reasons, the sample selected in the present paper is far from perfect but likely to contain 
unphysical pairs which would lead to wrong estimates of M/L. Therefore, the exclusion of 
true pairs and the contamination of optical pairs must be taken into consideration for M/ L 
determination. 

Fortunately, the possibility of true-pair exclusion is independent of i?p or Vp, as the 
selection criteria do not depend on them. Hence, the probability distribution of pairs in 
Rp-Vp space, which is to be compared with the observed pairs, can be expressed as 



where / is a constant corresponding to the fraction of true pairs out of observed pairs. 
Clearly the first term on the right side expresses the contribution of true binaries, and the 
second term describes the contamination from optical pairs. Probabilities p, Pbin, and Popt 
are normalized so that 



where the integrations are performed from to 400 kpc for Rp, and from to 400 km s 



The possibility of mis-identification of optical pairs is proportional to the number 
density of galaxies. It is generally known that the distribution of galaxies in the Universe 
is not uniform but shows strong clustering, which is usually described in terms of the 



p[Rp, Vp\M/L, f] = /pbin[it:p, Vp\{M/L)] + (1 - f)popt[Rp, Vp], 



(13) 




(14) 



for Vp. 
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two-point correlation function (e.g. Peebles 1993). With this function the probability 
distribution of optical pairs in the Rp-Vp phase space can be written as 

Popt[i?p, ^p] dR^dV^ (X [1 + e(r)] i?p dR^dVp, (15) 

where ^(r) is the two-point correlation function, and this is usually written in the form of 

= (16) 

The two-point correlation function is determined well in the scale of 10 Mpc but less certain 
in the scale of 1 Mpc. Hence in the following analysis we consider two cases, no clustering 
case with g = and clustering case with q — 1.8 and Tq = 10 Mpc (Peebles 1993), and see 
how the M/L estimates depend on the clustering effect. The separation r can be calculated 
from the projected separation and the velocity difference by assuming the Hubble constant 
of 50 km s"^ Mpc"^. Note that in any case the probability distribution of optical pairs in 
the -Rp-V^p space is independent of the M/ L ratio of galaxies. 



3.3. M/L Determination 

To evaluate the mass-to-light ratio and the fraction of true pairs we make use of the 
maximum- likelihood method for M/L and /. The probability for finding a pair at Rp and Vp 

in the projected phase space is proportional to p[Rp. Vp\M/L, /], and hence the logarithmic 
likelihood for finding all observed pairs at their observed positions in the projected phase 
space is expressed as the summation of the probability for finding each pair at its position. 
Therefore, the logarithmic likelihood of M/L for observed pairs can be written as 

log£(M/L, /) = ^n{Rp, Vp) logp [Rp, Vp\M/L, /], (17) 

where n{Rp, Vp) denotes the observed number of pairs having Rp and Vp, and the 
summation is done over the whole projected phase space {Rp less than 400 kpc and Vp less 
than 400 km/s). Evidently n{Rp,Vp) is integral as long as the values of Rp and Vp are 
determined with sufficient accuracy. However, for the pairs we consider here Rp and Vp 
have uncertainties, and the uncertainty in Vp is particularly crucial for M/L determination 
because the M/L estimator depends on V^. Therefore, we treated each observed pair as a 
Gaussian distribution spread in the direction of Vp, and then n{Rp, Vp) is given as 

n{Rp, Vp) dRpdVp (xY,9i dRpdVp, (18) 

i 

where 

^. = (2vra0-^/^exp(^^^^^). (19) 
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Here and (7j denote the observed Vp, and the uncertainty for Vp for ith pair, respectively. 
Note that n is normahzed so that J J n dRpdVp — Ni^t, where A^tot is the total number of 
pairs. 

Since the probability distribution of physical pair in the Rp-Vp phase space, 
Phin[Rp:Vp\{M/L)], cauuot be expressed analytically due to the projection effect, we 
performed Monte-Carlo simulations to evaluate pbin- The distribution of one million pairs 

in the projected phase space were simulated assuming random orientation of orbital planes 
with respect to the line of sight, and then the logarithmic likelihood (equation [16]) was 
calculated in the parameter space of / and M/L. 

Figure 2 shows the likelihood contours for sample 1 (57 pairs) for the case of q = (no 
clustering for optical pairs). The thick lines arc for /? = (isotropic orbit) and dotted lines 
are for j3 = —oo (circular orbit). The figures show that the M/L estimates are not affected 
strongly by the assumed orbital parameters. The best estimates of M/L are 35I5 for /3 = 
and 28l| for /3 = —00, with true binary fraction / of O.SSlg;?^ for both cases (the error bars 
denote the 68 % confidence level) . The results for the true binary fraction / indicates that 
most of pairs in sample I are likely to be bound. The expected number of optical pair is 
about 7 out of 57 pairs, which is comparable to the number of pairs which appear in Figure 
lb above the envelope corresponding to {M/L)p of 20. 

On the other hand, figure 3 shows the hkelihood contours for sample I hke figure 2, 
but for q — 1.8 (clustering for optical pairs). The contours for two orbital models are 
shown, and again one can see that the weak dependence of M/ L on the orbital parameter. 

However, the true pair fractions / arc quite different from those for no clustering cases, as 
we obtained / = 0.71l'^o;i5 {P = 0) and / = 0.73to;l5 {P = —00) for clustering case. This is 
because the expected number of optical pairs with small Vp is much larger than that for no 
clustering case, and hence more number of galaxies with small Vp arc regarded as optical 
pairs. However, the best estimates of M/L are 36I4 for /3 = and 30^3 for P = —00, which 
are fairly close to those for no clustering cases. 

In figure 2 and 3 the best estimates of M/L change little depending on /. This is 
because the M/L is essentially determined by the pairs which lie below the envelope in 
figure lb: in the most range of / (e.g. / less than 0.95) the pairs far above the envelope are 
hkely to be optical pairs, and have little effect on the M/L determination. On the other 
hand, if / is set to be almost unity, the most likely M/L could become as high as 100 
to explain the pairs with high (M/L)p without optical pairs. However, the hkelihood for 
finding / of almost unity is very small compared to that for / between 0.6 to 0.9, for in that 
case the concentration of pairs to small {M/L)p seen figure 1 cannot be explained at all. 
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In order to test whether M/L and / obtained above can reproduce the distribution 
of observed pairs, we simulated distribution of model pairs in the R^-V^ phase space using 
the best-fit parameters. Figure 4 shows the simulated distribution of 57 modeled pairs 
with / = 0.88 and M/L = 35 assuming isotropic orbits for true pairs and no clustering 
for optical pairs. The figure resembles well the observed distribution in figure lb in many 
aspects; small fraction of pairs with extremely high (M/L)p, concentration of pairs to small 
Vp, and the envelope corresponding to (M/L)p ~ 20. This simulation confirms the vahdity 
of the results. 

In order to see if these results depend strongly on the sample we used, we also 
performed the same analysis for sample II (109 pairs), which are selected with weaker 
isolation criterion (a = 1.5). As seen in Figure la, this sample is likely to contain more 
number of optical pairs than sample I due to the weak selection criterion. In fact, the 
resultant value of / is O.SOloiog clustering case {q — 0) and 0.62lo;ii for clustering 

case (g = 1.8), with P = for both cases. However, the best estimates M/L are 35^3 
{(3 = 0) and 28li (/? = -oo) for no clustering cases, and 36^^^ {(3 = 0) and 30t| (/? = -oo) 
for clustering cases. These results remarkably agree with those for sample I, indicating that 
the results does not depend strongly on the samples. The results for sample II as well as 
those for sample I are listed in Table 2. 

However, we would like to note that the results might be changed if we take the other 
limit of /9; /9 = 1 corresponding to radial orbits. In this case the best M/L for sample I 
was found to be 421^^ (g = 0), and so the M/L would exceed 50 within 68 % confidence 
level. Yet the assumption of /? = 1 seems too radical, because in this case the pairs suffer 
from direct encounters that will probably lead them to mergers. In order for bound pairs to 
survive for many orbital periods, their orbits must be elliptical at least to some degree, and 
so P cannot be too close to unity. On the other hand, perfectly circular orbits (/3 = — oo), 
are also unlikely because this requires fine tuning of orbital parameters. Therefore, the 
results with (3 — presumably represent best the mass-to-light ratio of true pairs. 



3.4. M/L for pure spiral pairs 

In the M/ L determination above, we did not consider the variation of M/ L among the 
sample galaxies. The M/ L of galaxies are, however, usually considered to vary depending 
on the type of galaxies. Indeed, previous binary galaxy studies claimed larger M/L for 
ellipticals than that of spirals (e.g., Schweizer 1987). Therefore, it is interesting to study 
the M/ L ratio of galaxies for a specific type. 
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The 57 pairs in sample I consist of spirals, SOs, ellipticals and some others such as 
peculiars. The dominant type among them is spirals, which occTipics a fraction of 70% in 
sample I, and hence we try to estimate M/ L for spiral galaxies. In particular we concentrate 
on 'pure' spiral pairs which consists of two spiral galaxies later than Sa, because a pair 
of a spiral and an SO or elliptical do not necessarily reflect the M/L of spiral galaxies. 
We selected 30 pure spiral pairs out of 57 pairs in sample I (hereafter sample III), and 
performed the same analysis described above. The likehhood contours for no clustering case 
{q — 0) and for clustering case {q — 1.8) are plotted in figure 5 and 6, respectively. Most 
likely value for no clustering case is found to be 15l3 for /3 = 0, and 12jl| for /3 = — oo, 
which are compared to the results for the 57 pairs with mixed types, 35 for /3 = or 28 
for /3 = — oo. The best estimates of M/L for clustering case are similar to those for no 
clustering case whereas the true pair fraction / is relatively smaller (see tabic 2). In both 
cases the difference in M/L between sample I and III is significant, being above the 3o" 
level. Therefore, we conclude the difference is real, and that M/L for spirals are smaller to 
ellipticals or SOs. This is consistent with previous studies of binary galaxies, although the 
M/L obtained here are somewhat smaller than those from previous studies. 



4. Discussion 

4.1. M/L Dependence on Separation 

In the previous sections, pairs are approximated as two point masses. However, real 
galaxies have finite size, and it could be as much as 100 kpc if dark halos extend well beyond 
the optical disks. In this case, the approximation of point masses is not vahd, and M/L 
obtained above could be underestimation, particularly for pairs with small separations. 
Here we investigate the dependency of M/L on the separation of pairs, and test if the 
assumption of point masses is reasonable for the present samples. 

We divided 57 pairs in sample I into 3 subgroups depending on the separation. The 
three subgroups consist of 27 pairs with < -Rp < 100 kpc, 12 with 100 < i?p < 200 
kpc, and 18 with 200 < kpc, respectively. The M/L ratios for three subsample were 
obtained in the same manner described in the previous section. Assuming /3 = and 
g = 0, we obtained the best estimates of M/L with la errors to be 36l5°, 37^^^, and 25^i3, 
respectively. The error bars are increased compared to the results in Section 3 because the 
number of galaxies in a sample is reduced, figure 7 shows the M/L dependence on the 
mean separations of pairs. The figure demonstrates that the M/L ratio is almost constant, 
and that the variations are within the la error bars. If the density distribution of dark halo 
is proportional to at a large radius, the M/ L increases linearly with radius, and if it is 
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proportional to as suggested by recent simulations (Navarro et al.l996), M/L increases 
with radius logarithmically. However, figure 7 shows no tendency of increasing M/L at 
large radii. Therefore, the halos of galaxies as luminous as the Milky Way Galaxies may be 
truncated within 100 kpc. The indication of constant M/L beyond 100 kpc is consistent 
with previous studies of binary galaxies. According to Peterson (1979b), the M/L ratio is 
gradually increasing with radius aX R < 100 kpc, but remains constant beyond it. Schweizer 
(1987) also showed that M/L does not increase beyond 100 kpc. 

The M/L estimate for spiral galaxies is also consistent with previous studies. Schweizer 
(1987) obtained M/L of 21 ±5 (V band, absorption corrected) with the sample whose mean 
separation is about 90 kpc. Peterson (1979b) obtained spiral galaxies' M/L of 35 ± 13 
[Hq ~ 50 km s~^ Mpc~^) based on 39 pairs with their mean separation of 110 kpc, and 
Turner (1976b) also obtained M/L for spiral galaxies of ~ 35. Note that in the 70's the 
correction for the galactic and internal absorption were not usually made, and this partly 
explains smaller M/L in the present paper. The mean of absorption correction in our 
sample galaxies is about 0.4 mag, which reduces M/L by ~ 30 %. Therefore, if similar 
amount of absorption correction were made, the studies by Peterson (1979b) or Turner 
(1976b) would give M/L of 23~27, which are close to our results of M/L. Although the 
M/L obtained in the present paper is not significantly different from previous studies, we 
would like to emphasize that the mean of absolute separations (not the luminosity corrected 
separation R^) of sample III in this paper is ~ 206 kpc, which is almost twice of those in 
previous studies. Nevertheless, the M/ L ratio does not show any tendency of increase with 
increasing the separation when compared with previous studies. 



4.2. Dcirk Halo Extent of Spiral Galaxies 

Since the mass of spiral galaxies' optical disk can be estimated from rotation curves, 
we can compare the total M/L with optical disk M/L, and discuss the extent of dark 
halos of spiral galaxies. We define (M/L) ^25 as the ratio of the enclosed mass within i?25 
to the total luminosity, where R25 is a radius at which the surface brightness becomes 25 
mag per arcsec^. The enclosed mass can be estimated from the HI velocity line width VFhi- 
Assuming that VFhi corresponds to twice of the rotation velocity, we obtain 

(M/L)h.. = (20) 

The central surface luminosity of spiral galaxies is constant, about 22 mag per arcsec^ 
(Freeman 1970). If this applies to the spiral galaxies in the binary sample, R2Z corresponds 
to about 3 times disk scale length d, and (M/L) ^^25 roughly approximates the mass-to-right 
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ratio of the disk. We calculated (M/L) ^^25 of spiral galaxies in sample I for which R25 
and Whi are available. R25 and Wm were taken from de Vaucouleurs et al.(1991), and 
Huchtmeier and Richter (1989), respectively. The values of {M/L)ji25 range from 2 to 18 
with the average of 7. Note that this M/L is consistent with the previous studies for disk 
M/L; for example, Faber & Gallagher (1979) obtained the M/L of about 5 {Hq = 50 km 
s~^ Mpc^^). This M/L is compared to the total M/L obtained in the section 3, M/L of 12 
~ 16. The total M/L is somewhat larger than the disk M/L, and the difference is almost 
3(7 level (see figures 5 and 6). This difference is of course due to the dark halo, and this 
indicates that under the maximum disk assumption the contributions of dark halo and 
the optical disk to the total mass of galaxies are comparable. However, the assumption 
of maximum disk is still controversy. If the disk mass is smaller than that indicated by 
the maximum rotation velocity within the optical disk, the dark halo could dominantly 
contribute to the total mass. 

We can also estimate the extent of dark halos by comparing the total M/ L with disk 
M/L. If a flat rotation curve is assumed, the value of M/L increases linearly with radius. 
In this case, the resultant M/L of 15 implies that the typical halo extends 15/7 ~ 2 times 
i?25- If we adopt disk M/L of 5 according to Faber & Gallagher (1979), the halo extent 
is 15/5 ~ 3 times i?25- Therefore, the typical dark halo size maybe about 6 to 9 times 
disk scale length, if R25 is ~ 3 times disk scale length d. This is to be compared with the 
size of optical disk, which is about 4 ~ 5 times disk scale length (van der Kruit & Searle 
1981). Hence, if the rotation curve is perfectly flat out to the radius at which the halo 
mass distribution is truncated, the halo size may be about 2 times larger than optical disks. 
This is, of course, not exactly true when the rotation curve is not completely flat, which is 
preferred by the recent simulations. In this case, the halo size may be somewhat lager, but 
it cannot exceed a few hundred kpc as indicated by the M/L constancy. 

The halo size indicated here is somewhat smaller than those in previous studies, but 
quite consistent with recent investigations. For instance, a number of declining rotation 
curves, which may be fitted even with a Keplerian, are found recently (e.g., Jorsater & 
van Moorscl 1995; Oiling 1996; Honma & Sofue 1996). Honma & Sofuc (1997) showed 
that such declining rotation curves arc not uncommon by considering the observational 
uncertainty. These rotation curves studies arc generally based on the observation of HI, 
which are usually observed out to 5 ~ 10 times disk-scale length. Therefore, the fact that 
the declining part of rotation curves were found is consistent with the present results. 

We are grateful to Y. Sofue for his supervision, and to Y. Tutui and J. Koda for fruitful 
discussion. This work was financially supported by the Japan Society for the Promotion of 
Science. 
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figure captions 

Figure 1. Distribution of selected pairs in the Rp-Vp phase space. Figure la is for the pairs 
selected with a = 1.5, and lb is for the pairs selected with a = 2.5. The dotted curve in 
figure lb corresponds to the constant {M/L)p of 20. 

Figure 2. Likelihood contours in the parameter space of M/L and / for sample I. As for 
the optical pair distribution no clustering is assumed (g = 0). Solid lines are for f3 = 
(isotropic orbits), and dotted lines are for (3 = — oo (circular orbits). Three contours for 
each case correspond to 68, 95 and 99% level, and the crosses denote the peak of the 
likelihood. 

Figure 3. Likelihood contours same to figure 2, but for q = 1.8 (clustering case for optical 
pairs). 

Figure 4. Simulated distribution of 57 pairs which are to be compared with figure lb. The 
distribution is calculated with the best fit parameters obtained in the Section 3. The 
dotted curve corresponds to (M/L)p of 20. 

Figure 5. Likehhood contours same to figure 2 (no clustering case), but for sample III (30 
pure spiral pairs). 

Figure 6. Likehhood contours same to figure 3 (clustering case), but for sample III (30 pure 
spiral pairs). 

Figure 7. M/L for three subgroups against the mean separation (see text for the sample). 
The error bar for R denotes the la deviation of pairs in the sample. 
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